Neutral pion threshold production at Q 2 = 0.05GeV 2 /c 2 
and Chiral Perturbation Theory 
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New data are presented on the p(e, e'p) 71 " reaction at threshold at a four-momentum transfer of 
Q 2 = 0.05 GeV 2 /c 2 . The data were taken with the three-spectrometer setup of the Al Collaboration 
at the Mainz Microtron MAMI. The complete center of mass solid angle was covered up to a center 
of mass energy of 4MeV above threshold. Combined with measurements at three different values 
of the virtual photon polarization e, the structure functions cr-r, gl, o~tt, and otl are determined. 
The results are compared with calculations in Heavy Baryon Chiral Perturbation Theory and with a 
phenomenological model. The measured cross section is significantly smaller than both predictions. 
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Introduction. Threshold electromagnetic pion pro- 
duction is a fundamental process since the pion is a 
"Goldstone Boson", reflecting the spontaneously broken 
chiral symmetry of QCD[|l). In the chiral limit, where 
the quark and pion masses vanish, the s wave produc- 
tion amplitudes of neutral pions vanish. However, the 
explicit chiral symmetry breaking due to the small but 
finite quark mass (m u « 5MeV/c 2 , w 9MeV/c 2 ) and 
finite pion mass render these amplitudes finite. 

Calculations of these observables are performed by an 
effective field theory called Chiral Perturbation Theory 
(ChPT) p|] , which is generally in good agreement with 
experiment 0. The systematic application of ChPT to 
reactions involving heavy baryons by Ref. || (Heavy 
Baryon Chiral Perturbation Theory, HBChPT) has been 
generally successful in describing 7r — N scattering and 
electromagnetic pion production from the nucleon. 

In recent years there has been a considerable experi- 
mental effort to test this theoretical approach^]. Of spe- 
cific interest to this work, photo production experiments 
were performed at Mainzg | § and at SAL@, § and 
showed an impressive agreement with the predictions of 
Ref. ||. These experiments were extended to finite pho- 
ton four- momentum transfer — Q 2 via electro production 
at NIKHEF§, [kJ and MAMI@ at Q 2 = 0.1 GeV 2 /c 2 , 
which is believed to be the limit of the predictive power 
of HBChPT. Nevertheless, the results were in reasonable 
agreement with the calculations [fl2||. 

In this paper we present a measurement at a value of 
Q 2 = 0.05GeV 2 /c 2 , half way between the photon point 



and the existing electro production data. The present 
data cover the complete azimuthal angle <f> and allow 
for a Rosenbluth separation to extract the unpolarized 
structure functions ax, and otl- err was observed 
to be so small that only upper limits were obtained. 
Since the former experiments agree with the predictions 
of HBChPT, the clear disagreement between our results 
and these calculations are surprising. 

Formalism and Kinematics. In the one-photon ex- 
change approximation, the unpolarized electro produc- 
tion cross section can be written as (see e.g. Ref. 
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with the virtual photon flux T, photon energy lo* , and 
transverse polarization e, er, — eQ 2 /u>* 2 . The pion emis- 
sion angles 9* and </>* are relative to the momentum 
transfer q and the scattering plane. Variables in the 
photon-proton center of mass frame are denoted with an 
asterisk. 

In the threshold region, where only s and p waves have 
to be taken into account, the cross section can be further 
decomposed as 

o- T (6*) = (pi / k;) (a + b cos e: + c cos 2 e;), 

o- L {6%) = (p*Jk;)(A' + B' cos 91 + C" cos 2 91), 
vtl(K) = (p*/K)(Dsin9* w + Esm6*cos6l) , 
cr TT {0*) = (p:/k*)(Fsm 2 9:), (2) 
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where p^/k* is the ratio of pion CM momentum and 
photon CM equivalent momentum. The angular coeffi- 
cients are given by two s wave and five p wave multipole 
combinations fll2| 
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At threshold, the s wave multipoles Eq + and Lq + are 
real and their energy dependence is governed by the uni- 
tary cusp caused by the two step 7*p — > ir + n — > ir°p 
amplitude |l4|. The p wave amplitudes P, are real and 
proportional to the pion center of mass momentum p* . 
The coefficient F is too small to be extracted in this 
experiment, so only the p wave combination P 2 3 = 
(|P 2 | 2 + |P 3 | 2 )/2 can be determined. 

Experiment. The experiment was performed at the 
three-spectrometer setup of the Al Collaboration at 
MAMI (see Ref. |15 for a detailed description). A liq- 
uid Hydrogen target with a length of 49 mm and target 
walls of 10 /im Havar was used with beam currents of 
up to 30 /iA, i.e. at a luminosity of 39MHz//ib. For the 
detection of the scattered electron, spectrometer B with 
5.6 msr solid angle at 15% momentum acceptance was 
used for the two forward angles, while spectrometer C 
with 22.5 msr solid angle at 25% momentum acceptance 
was used for the backward angle. 

For the detection of the recoil proton, spectrometer A 
with 21 msr solid angle and a momentum acceptance of 
20% was used for all settings. Table | shows the central 
angle and momentum settings of the spectrometers. Due 
to the kinematical focusing by the Lorentz boost, the 
full center of mass solid angle was covered within each 
setting up to a center of mass energy of AW = 4MeV 
above threshold. 

All spectrometers were equipped with a detector pack- 
age consisting of four layers of vertical drift chambers 
for position and angle reconstruction, and two layers of 
scintillator detectors for time-of-flight measurement and 
trigger. For the electron arm, the trigger was defined by 
the coincidence of both scintillator layers, while the low 
energy protons were already stopped in the first scintilla- 
tor layer of spectrometer A. A halocarbon gas Cercnkov 
detector with a detection efficiency of 99% in this energy 
range was used for n~ je~ separation, but not as part of 



TABLE I: Kinematical settings for the central trajectories 
of the spectrometers. The four-momentum transfer for all 
settings is Q 2 = 0.05GeV 2 /c 2 
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the online trigger decision. The electron detection effi- 
ciency and calibration was checked by a measurement of 
elastic electron scattering. 

Analysis and Error Estimate. The coincidence time 
between electron and proton spectrometer, corrected for 
the flight path in the spectrometers, was determined with 
a resolution of 1.8 ns FWHM. This was limited by the un- 
certainty in the flight path due to multiple scattering of 
the low energy protons in the detector package of spec- 
trometer A. 

After a cut on the coincidence time, the missing 
mass was determined by the four-momentum subtrac- 
tion of incoming and outgoing particles, i.e. by rn 2 miss — 
{e-in +Pin - &out - Pout) 2 ■ A missing mass resolution of 
2.0MeV/c 2 was achieved, again limited by the multiple 
scattering of the low energy protons. 

The phase space was calculated in a Monte- Carlo sim- 
ulation including resolution, multiple scattering, energy 
loss, and radiative corrections, the later based on the for- 
mulas of Ref. @. 

The background contribution was determined by a 
cut in the coincidence time spectrum. Since the mod- 
est missing mass resolution allowed only a loose cut of 
±3 MeV/c 2 around the 7T° mass, the signal to background 
ratio was « 1 averaged over the acceptance. The back- 
ground, however, could be determined with a timing cut 
of 10 times the width of the cut for the true events, lead- 
ing to no significant contribution to the statistical error 
by the background subtraction. 

The systematic error is dominated by different effects 
at low and high energies. For the lowest energy bin, 
the calibration of the momentum of the electron arm is 
crucial, since assuming a too high electron momentum 
leads to a smaller photon momentum and a smaller total 
CM energy, so that some of the events are shifted below 
threshold. This error was estimated by varying the cen- 
tral momentum of the electron arm in the analysis. For 
the lowest energy bins we assumed Sa/a = 20% error at 
AW = 0.5 MeV and 10% at AW = 1.5 MeV. 

For the higher energy bins, the acceptance of the pro- 
ton arm is not uniform, and a few angular bins are mul- 
tiplied by large phase space correction factors. This con- 
tribution to the systematic error can be estimated by 
varying software cuts on the acceptance. 

Compared to these two effects, the contribution of the 
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FIG. 1: Differential cross sections for the first 4 MeV above 
threshold for the virtual photon polarization e = 0.72. The 
solid line represents a fit with the assumption of only s and 
p waves contributing, the dashed and dash-dotted lines rep- 
resent the predictions of HBChPTlll and MAID ill]. 



efficiency and other corrections can be neglected. Over 
all, we estimated an error of da /a = 5% (3%) at AW — 
2.5 MeV (3.5 MeV). 

After determining the two- fold cross section a (8* , </>*) 
for the three epsilon points, a fit with the known </>* de- 
pendence (Eqn. |l|) was performed in each 8* bin to ex- 
tract a Q (8*) = a T (6*) + e L a L (8*) and a TL {8l). 

Results. Fig. [l] shows the differential cross section for 
four energy bins at the medium epsilon point in com- 
parison with the predictions of HBChPT|l2| and a phe- 
nomenological Model (MAID) ]l7[) . The complete data 
set can be obtained from our web site. Fig. || shows 
the transverse-longitudinal interference structure func- 
tion extracted as weighted average from all three settings. 

In order to extract the multipole amplitudes, a fit us- 
ing Eqn. |§| with the assumption of constant real s wave 
multipoles (i.e. Eq+, Lq + constant in energy) and real p 
wave multipoles proportional to the pion CM momentum 
(Pi = p*Pi with Pi constant in energy) was performed, 
indicated by a solid line in the figures. These assump- 
tions are not exactly valid, since a variation of w 5% in 
the total s wave amplitude is expected due to the unitary 
cusp|l^, but the statistical significance of the data is 
not sufficient to resolve this variation. The extracted fit 
parameters are presented in the first row of table [fi|. 

It is important to note, that a least squares fit cannot 
give a real picture of the sensitivity of the data to the 
multipole amplitudes. The data set is dominated by the 
systematic error and thus violates a fundamental precon- 



FIG. 2: The transverse-longitudinal interference structure 
function, determined as weighted average of all three settings 
for e. Assignment of lines as in Fig. pL 
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FIG. 3: The total cross section a to t versus Q 2 , at a value of 
e — 0.8. The solid (dashed) line is the prediction of ChPT 
(MAID), data points at Q 2 = and 0.1 GeV 2 /c 2 from |, 0. 



dition of Gaussian distributed errors for a % 2 fit. 

To compare the present result with former experiments 
in a consistent way, the same fit was performed to the 
data sets at Q 2 = 0.1 GeV 2 /c 2 , which are believed to be 
more or less consistent with HBChPT on the cross sec- 
tion level. The fit results are included in table [n] and 
differ from the multipole amplitudes quoted in the cor- 
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TABLE II: Extracted multipole amplitudes in comparison 
with the threshold values of HBChPT[jt2| and MAID M . The 
AmPS[[lo| value for [Lo+] was extracted from their value for 
do ~ £l\Lo+ | 2 . For the AmPS § fit Lo+ was fixed, since no 
Rosenbluth separation was performed. 
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responding references^, O], which were extracted using 
model assumptions. 

Fig. shows the unpolarized total cross section otot — 
J (ctt + clCl)^ as a function of Q 2 to further illustrate 
the discrepancy between the currently available data and 
the calculations in a model independent manner. While 
the photo production data could be described by both, 
HBChPT and MAID, the strongest deviation appears at 
the Q 2 value of the present experiment, while already at 
Q 2 = 0.1GeV 2 /c 2 , where the parameters of HBChPT 
were determined, a deviation from the HBChPT calcu- 
lation of up to 20% appears. The disagreement with the 
MAID model is even larger. 

From the fit parameters quoted in table |H| one can see, 
that the major part of the discrepancy comes from the 
large p wave multipole combination P 2 3 . This is a serious 
problem for the prediction of HBChPT, since this combi- 
nation is already fixed by the photoproduction data and 
cannot be adjusted by varying the free parameters of the 
calculation to the current order. The strong curvature 
in the Q 2 dependence of Fig. || clearly indicates, that 
higher orders in Q 2 have to be included in the p wave 
description. 

In order to check the surprising result of our experi- 
ment, we repeated our differential cross section measure- 
ment for the highest epsilon (0.92) in an independent 
experiment. The result of this experiment agreed with 



the present experiment within the errors. In addition, 
since these results rely on the quality of the simulation 
and phase space integration program, a new, completely 
independent code was written to check the phase space 
integration. 

In summary, it appears that there is a significant dis- 
crepancy between HBChPT and phenomenological mod- 
els on the one hand and the experimental data on the 
other hand. However, an inconsistency in the different 
electro production data sets, which were all taken in sep- 
arate experiments, can not be excluded. As mentioned 
above, a repeated measurement at Q 2 = 0.05GeV 2 /c 2 
and e = 0.92 at MAMI confirmed our result. We plan 
to further explore the Q 2 dependence in a future experi- 
ment. An independent measurement at the Jefferson Lab 
in the same Q 2 range is also planned Jl8|. 
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